Quantum Tunneling and Quantum-Classical Transitions in Large Spin Systems 
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We have studied quantum tunneling of large spins in a biaxial spin system and in a single-axial 
spin system with a transverse magnetic field. The asymptotically exact eigenvalues and eigenstates 
of the spin systems are obtained by solving the Mathieu equation. When the height 4q of energy 
barrier formed by the anisotropy and the magnetic field exceeds some critical values 4qiT{r), the 
ground states and the excited states split due to tunneling of large spins. We also have presented 
the phase diagram of these spin systems in the extremely weak spin-phonon coupling limit. The 
crossover from thermal to quantum regime is second order. With further decreasing temperature, the 
first order phase transition occurs from a thermally-assisted resonant tunneling to thermal regime 
as q < gir(s). Our theory agrees with recent experimental observations well. 
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In recent years quantum tunneling of magnetization 
(QTM) in small magnetic particles has been widely in- 
vestigated in both theory and experiment because of its 
fundamental interest in exploring the transition between 
classical and quantum physics [1]. It gives evidence of 
quantum mechanical behavior in macroscopic systems. 
Such phenomena have been observed experimentally in 
ferritin [2], the molecular nanomagnets Mni2-ac [1, 3, 4], 
Fes [5, 6], BaFeCoTiO [7], etc at very low temperatures. 
Theoretical calculations of such effects were performed 
for the ferromagnets [1, 8] and antifcrromagnets without 
and with the excess spin [1, 9] by using path integral 
and WKB methods, etc.. However, these approaches de- 
pend on the semiclassical treatment and did not present 
the energy spectrum of the spin systems, which is nec- 
essary to give the complete solutions of spin tunneling 
problems. Especially, the quantization of spin levels be- 
comes very important to explain well the experiments of 
quantum tunneling at very low temperatures [4] . In very 
recent Rcfs. [10-12], the spin systems arouse new interest 
again because they provide examples exhibitting first or- 
der phase transitions, which were not known previously. 
The energy spectrum of the spin systems can also help 
us to understand the origin of the first order phase tran- 
sition. In this paper, we investigate a biaxial spin system 
without the magnetic field and a single-axial spin sys- 
tem with a transverse magnetic field which are generic 
for spin tunneling problems studied by different methods 
[1, 8, 10-12] in the framework of the Schrbdinger's pic- 
ture of quantum mechanics, so that the energy spectrum 
as well as the physical properties of the spin systems are 
obtained. 

We first consider QTM in the biaxial spin system with- 
out the magnetic field described by the Hamiltonian [1, 



11] 



H — -AS 2 — BS 2 , 



(1) 



where the anisotropy constants A and B are taken 
to be positive. So the ground state of the system 
corresponds to spin S pointing the positive or neg- 
ative x axis. To diagonalize the Hamiltonian (1), 
we first write its matrix representation in the basis 
\s,m >, which has the properties: S 2 |s,to >= s(s + 
l)|s, m >, S z \s, to >= m|s, to > and S±\s, to >= (S x ± 
iS y )\s,m >— -\/s(s + 1) — m(m± l)|s, to ± 1 >, where 
s is the spin quantum number to be taken as an integer 
[13], m = —s, — s + 1, • • • , s and the unit of h = 1 is used. 
Let E and ^ m be the eigenenergies and the eigenstates 
of H , respectively, then we have 



(to 



+ y/[s(s + 1) - (to + l) 2 ] 2 - (to + l)^ m+2 (2) 
+4[f + ±s(s + 1) - (A + i)rn 2 ]* m = 0, 

where A = -p. For a large spin system, s is a sufficiently 
large number, i.e. s ^ 1. Define $ m = (— l)^]^ r m , 
where [^] denotes the integer part of ^^then Eq. (2) 
becomes 

, N d 2 $ d$ r E 1 , . ,,11, 



in the large- s limit. Here x 



(3) 

and only the 



leading terms are remained. Taking the transformations 
<f> = (1 — x 2 )~*y(x) and x = sint and substituting them 
into Eq. (3), we finally obtain the well-known Mathieu 
equation 



dt 2 



+ [A( 9 )-2 g cos(2i)]t/ = 0, 



(4) 
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which describes the motion of a particle with the mass 
5 in the cosine potential. Here the characteristic values 
A(q) = —jjj + 2q and q = \Xs(s + 1). Obviously, when 
q = 0, Ai TO i(0) = m 2 , then E m = —Bm 2 , which are 
the exact eigenenergies of (1) in the eigenstates \s, m > x . 
Now we have completed the mapping of the spin prob- 
lem onto a particle problem. The energy spectrum of the 
Hamiltonian (1) is determined completely by the Math- 
ieu equation. 

The Mathicu equation (4) has been studied in a large 
set of literature [14, 15] due to its physically basic im- 
portance. It is known that there exist periodic solutions 
of period nw, where n is any positive integer. However, 
the solutions relative to quantum tunneling problem dis- 
cussed here are only those even and odd ones with periods 
7r and 2ir, i.e., y =ce r (t, q) and se r (t, q), r = 0, 1, • • • , s. 
Assume that a r and b r are the characteristic values asso- 
ciated with the even and odd periodic solutions, respec- 
tively, then these characteristic values form a countable 
sequence, i.e. ao < b\ < a\ < b 2 < »2 < • • • < b s < a s 
for q > and a 2i (-q) = a 2i (q), b 2i (-q) = b 2 i(q), 
a2i+i(-q) = b 2 i+i(q) and b 2i+ i(-q) = a 2i+1 (q). We note 
that the approximately degenerate characteristic levels 
(a r ,b r ) are lifted remarkably as the energy barrier pa- 
rameter q exceeds the critical value qvrir). The higher 
the characteristic level is, the larger the critical value 
is, i.e. qvr{i) > 9it(?) f° r i > j- The value of qir(r) 
can be estimated by letting the height of the energy bar- 
rier be equal to the energy A r (0) of the particle in the 
vanishing energy barrier, i.e. 4qir(r) = r 2 , which is a 
good approximation for large r. This shows clearly that 
when q > gir(r), the splitting of the characteristic lev- 
els (a r ,b r ) is induced by tunneling of the particle. As 

q is larger than q 2T {r) w 4(^+1)2 I 16 ]' the levels br 
and a r _i degenerate approximately due to very high en- 
ergy barrier. Therefore, in order to observe experimen- 
tally tunneling splitting of level r, q must locate between 
qirir) and q%r(r). We also note that tunneling of the 
particle mainly occurs in those characteristic levels rang- 
ing from the highest level to one with characteristic value 
a s — s 2 +0(s) as q < qir{s), which is important to evaluat- 
ing approximately the critical temperature from thermal 
to quantum regime below. 

For the Hamiltonian (1), its lower excited states corre- 
spond to the higher characteristic levels of the Mathieu 
equation. In order to see the tunneling splitting of the 
degenerate ground states, q must exceed the critical value 
<?it(s) ~ \s 2 , i- e - Ait(s) ~ 1 — -, which coincides with 
the result reported previously [11] when s is very large. 
When q > q%r(m), we obtain tunneling splitting of the 
ground states as well as the excited states [14] 

AE m = B(a m - b m ) = B(2y/Xs(s + 1) - m) + 0(g - *). 

(5) 

When q — > 00, the ground state a s is singlet while all 
the excited states are twofold degenerate, which coincide 



with the degeneracy of the Hamiltonian (1) with B — > 0. 

Our another example is devoted to the single-axial spin 
system with a transverse magnetic field, which Hamilto- 
nian has form [1, 8, 10] 

H=-hS z -BS 2 x , (6) 

where B > and h > 0. The Hamiltonian (6) can be also 
diagonalized by the same method mentioned above. Its 
energy spectrum is determined by the Mathieu quation 
(4), too. However, in this case, x = cos(2t), A(q) = — ^ 

and q = 2/ V^+!) As h = 0, then A r (0) = r 2 and 
E m = —Bm 2 . So the ground states and the excited 
states of the Hamiltonian (6) correspond to the char- 
acteristic levels (a 2s ,b 2s ) and (a 2r ,b 2r )(r < s), respec- 
tively. According to that 4qiT(m) — (2m) 2 , we have the 
critical magnetic field hrrim) = — f m , m 3> 1. As 

q < hvrirn), there is no tunneling between the approxi- 
mately degenerate energy levels E m (i. e. a 2m and 62m) 
of the Hamiltonian (6). As h > /iit(s) = — , Bs , the 

ground states occur splitting, which is the same with that 
in Ref. [10] for very large s. We also note that h must be 
smaller than h c = 2Bs, at which the degenerate ground 
states coincide. As q » 1, we have 

AE m = \B(a 2m -b 2m ) = B(sj 2h ^ S + 1 ) -l m )+0( g -^). 

When q — > 00, all the energy levels a 2m and b 2m are sin- 
glet, which agree with the energy spectrum of the Hamil- 
tonian (6) with B — > 0. Because there is no Kramers' de- 
generacy in the Hamiltonian (6) , tunneling splitting also 
holds for half odd integer spin s [14]. 

Up to now, we have constructed the energy spectrum 
of the pure spin systems (1) and (6). However, the exper- 
iments of observing quantum tunneling were performed 
at low temperatures. So the influence of phonons or ther- 
mal activation on tunneling must be considered. Here, 
we assume the spin-phonon coupling is so weak that the 
energy spectrum of the spin systems is not changed too 
much. The transition between the ground state and the 
excited states can be completed by absorbing or emitting 
phonons. With the help of the Debye theory of a solid 
[17], we obtain approximately the critical temperature T c 
satifying 

3k B T c D(?f) w Bs 2 (8) 

as q < qir{s) and 

3k B T c D(f ) » B(a s - b r ) for (1) , , 
w \B{a 2s - b 2r ) for (6) ^ 

&sq 2T (r) <q< q 2 r(r+l) &ndq 2T (2r) <q< q 2T {2r+l), 
respectively. Here the Debye function D{x) w i7r 4 x 3 at 



2 



low temperature, fee and 9 are the Boltzmann constant 
and the Debye temperature, respectively. Obviously, the 
crossover from the high temperature (thermal) phase to 
the low temperature (quantum) phase is second order be- 
cause tunneling splitting of large spin continuously trans- 
forms into zero. As q < qiT(s), in order to see tunneling, 
large spin must absorb at least energy B(s 2 — Aq) and 
B(s 2 — q) for the models (1) and (6), respectively, pro- 
vided mainly by phonons or thermal activation. So we 
have the first order phase transition temperature T , 

3k B T D(^ ) « B(s 2 -4q) for (1). 

w B{s 2 -q) for (6). { ' 

The phase diagram of the spin systems (1) and (6) is 
shown in Fig. 1. 

Here we compare our theoretical results with recent ex- 
periments of quantum tunneling. The Hamiltonian (1) is 
believed to be a good description for the molecular nano- 
magnet Fes with spin s = 10 [5]. The parameters A — 
0.092K and B = 0.224K, then q = 11.295 < q 1T {s) = 25. 
Therefore, with decreasing temperature, the spin system 
enters the thermally-assisted resonant tunneling regime 
II from thermal regime I, as observed down to 0.067K. 
However, tunneling can not be seen possibly if the tem- 
perature is further lowered to below To, which would be 
of interest to be verfied experimentally. The Hamiltonian 
(6) has been found to be a good approximation for the 
molecular nanomagnct Mni2-ac with spin s = 10 [4]. For 
this nanomagnet, B = 0.68K, then /iir(s) = 3.24K. The 
parameter h = g±[iBH was changed from 6.38K to 8.93K 
in experiment, i.e. the magnetic field H was applied be- 
tween 5T and 7T. So the spin system fell into quantum 
regime III, where the two states of the fundamental dou- 
blet of the high-spin molecule are lifted. The transitions 
between the two states induced by a two-phonon process 
have been observed from 0.1K to 0.02K when H — 6.06T 
[4]. Obviously, as the spin systems just set in the thermal 
regime IV or q > <72t(s), tunneling does not be observed 
in experiments even at very low temperatures [18]. 

In summary, we have investigated quantum tunnel- 
ing of large spins based on simple models (1) and (6). 
The asymptotically exact spectrum of the spin systems 
is completely determined in the full range of the magnetic 
anisotropy and the magnetic field. The degeneracy of all 
the energy levels is removed due to QTM. The phase dia- 
gram presented here coincides with the experimental ob- 
servations. To our knowledge, the thermal regime IV in 
Fig. 2 was not predicted by previous theories. We think 
that the first order phase transition is not limited to the 
spin systems (1) and (6) and also exists generally in other 
spin systems, depending on whether the energy barrier 
and the ground states as well as the low-lying excited 
states cross. The method of solving the spin Hamilto- 
nians used above is also applied to the other large spin 
systems, such as the other symmetric ferromagnets, an- 
tiferromagnets and etc., which is in progress. 



This work is supported in part by grants from the Hong 
Kong Research Grants Council (RGC), the Hong Kong 
Baptist University Faculty Research Grant (FRG), the 
Sichuan Youth Science and Technology Foundation, the 
NSF of the Sichuan Educational Commission and the 
NSF of China (No. 19677101). 



[1] Quantum Tunneling of the Magnetization, edited by L. 
Gunther and B. Barbara, NATO ASI, Ser. E, Vol. 301 
(Kluwer, Dordrecht, 1995); Bo-zang Li and Wen-ding 
Zhong, in Aspects of Modern Magnetism (F. C. Pu et 
al eds., World Sci., Singapore, 1996)57. 

[2] D. D. Awschalom, J. Smyth, G. Grinstein, D. DiVin- 
cenzo, and D. Loss, Phys. Rev. Lett. 68, 3092 (1992); J. 

Tejada, X. X. Zhang, E. del Barco, J. M. Hernandez, and 

E. M. Chudnovsky, Phys. Rev. Lett. 79, 1754 (1997). 

[3] C. Paulsen, J. G. Park, B. Barara, R. Sessoli, and A. 

Caneschi, J. Magn. Magn. Mater. 140-144, 379 (1995); 

J. R. Friedman, M. P. Sarachik, J. Tejada, and R. Ziolo, 

Phys. Rev. Lett. 76, 3830 (1996); S. Hill, J. A. A. J. 

Perenboom, N. S. Dalai, T. Hathaway, T. Stalcup, and 

J. S. Brooks, Phys. Rev. Lett. 80, 2453 (1998). 
[4] G. Bellessa, N. Vernier, B. Barbara, and D. Gatteschi, 

Phys. Rev. Lett. 83, 416 (1999). 
[5] C. Sangregorio, T. Ohm, C. Paulsen, R. Sessoli, and D. 

Gatteschi, Phys. Rev. Lett. 78, 4645 (1997). 
[6] W. Wernsdorfer, T. Ohm, C. Sangregorio, R. Sessoli, D. 

Mailly, and C. Paulsen, Phys. Rev. Lett. 82, 3903 (1999). 
[7] W. Wernsdorfer, E. Bonet Orozco, K. Hasselbach, A. 

Benoit, D. Mailly, O. Kubo, H. Nakano, and B. Barara, 

Phys. Rev. Lett. 79, 4014 (1997). 
[8] J. L. Van Hemmen and A. Suto, Physica 141, 37 (1986); 

M. Enz and R. Schilling, J. Phys. C 19, 1765, L711 

(1986); E. M. Chudnovsky and L. Gunther, Phys. Rev. 

Lett. 60, 661 (1988); G.-H Kim and D. S. Hwang, Phys. 

Rev. B 55, 8918 (1997). 
[9] B. Barbara and E. M. Chudnovsky, Phys. Lett. A 145, 

205 (1990); I. V. Krive and O. B. Zaslavskii, J. Phys: 

Cond. Matter. 2, 9457 (1990); A. Chiolero and D. Loss, 

Phys. Rev. B 56, 738 (1997); S. E. Barnes, R. Ballou, B. 

Barbara, and J. Strelen, Phys. Rev. Lett. 79, 289 (1997). 
[10] E. M. Chudnovsky and D. A. Garanin, Phys. Rev. Lett. 

79, 4469 (1997). 
[11] J. -Q. Liang, H. J. W. Muller-Kirsten, D. K. Park, and 

F. Zimmerschied, Phys. Rev. Lett. 81, 216 (1998). 

[12] D. A. Garanin and E. M. Chudnovsky, Phys. Rev. B 56, 
11102 (1997); D. A. Garanin, X. Martinez Hidalgo and 
E. M. Chudnovsky, Phys. Rev. B 57, 13639 (1998); D. A. 
Garanin and E. M. Chudnovsky, Phys. Rev. B 59, 3671 
(1999); Gwang-Hee Kim, Phys. Rev. B 59, 11847 (1999); 
E. M. Chudnovsky, Phys. Rev. A 46, 8011 (1992). 

[13] For half odd integer spins, tunneling splitting of all the 
energy levels of the Hamiltonian (1) vanishes due to 
Kramers' degeneracy. 



3 



[14] M. Abramowitz and I. A. Stegun, Handbook of Mathe- 
matical Functions, (New York: Dover, 1965). 

[15] J. Meixner and F. W. Schafke, Mathieusche Funktionen 
und Sphdroidfunktionen, (Springer, Berlin, 1954). 

[16] The values of q2r{r) can be obtained by diagonalizing the 
Hamiltonian (1) in the basis \s, m > x , which also leads 
the Mathieu equation. qir{r) and q2r(r) coexist. 

[17] C. Kittel, Introduction to Solid State Physics, 7th ed. 
(Wiley, New York, 1996). 

[18] W. Wernsdorfer, E. Bonet Orozco, K. Hasselbach, A. 
Benoit, B. Barbara, N. Demoncy, A. Loiseau, H. pascard, 
and D. Mailly, Phys. Rev. Lett. 78, 1791 (1997). 



FIG. 1. The phase diagram of the spin systems (1) and (6): 
I and IV are the thermal regimes; II is the thermally-assisted 
resonant tunneling regime; III is the quantum tunneling 
regime. The transition from I to II and III is second order 
while the transition from II to IV is first order. 
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